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X 1 Introduction 

In this paper we consider the following family of Poisson brackets 

{M i5 Mj} = e ijk M fe , {M i; 7j} = e iifc 7*, {74, 7^-} = k e ijk M k . (1.1) 

Here Mj and 7« are components of 3-dimensional vectors M and T, e^k is the totally skew-symmetric 
tensor, k is a parameter. It is well-known that any linear Poisson bracket is defined by an appropriate 
Lie algebra. The cases k = 0, k > and k < correspond to the Lie algebras e(3), so(4) and so(3, 1). 
In this paper we assume n 7^ 0. 
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Bracket (jl.lj) has the two Casimir functions 

Ji = (M,r), j 2 = K |M| 2 + |r| 2 , 

where (-, •) stands for the standard dot product in M. 3 . Hence, for the Liouville integrability of the 
equations of motion only one additional integral functionally independent of the Hamiltonian and the 
Casimir functions is necessary. 

The simplest nontrivial class of quadratic homogeneous Hamiltonians of the form 

H = (M, AM) + (M, BT) + (I\ CT), (1.2) 

where A, B and C are constant 3 x 3-matrices, has many important applications in rigid body dynamics. 
We call Hamiltonians (jl.2|) having an additional polynomial integral integrable. We say that (jl.2)) 
possesses an additional integral of degree k if there is no non-trivial integral of degree less than k. In 
this paper all coefficients of both the Hamiltonian and the additional integral are supposed to be real 
constants. 

The case where the Hamiltonian (|1.2jl has a linear additional integral of motion has been investi- 
gated by Poincare [I]. 

There are two classical integrable cases, one found by Frahm-Schottky and one by Steklov, where 
the additional integral of motion is of second degree. It was proved in j2] that any Hamiltonian (Jl.2)) 
possessing an additional second degree integral is equivalent to one of these two cases. 

In 1986 Adler-van Moerbeke |Hj and Reyman-Semenov-tian Shansky ||] independently found a 
Hamiltonian of the form (jl.2)l with fourth degree additional integral. 

For the Frahm-Schottky, Steklov, and Adler-van Moerbeke-Reyman-Semenov-tian Shansky cases 
all matrices A, B and C are diagonal. This special subclass of " diagonal" Hamiltonians (jl.2j) was 
investigated by many authors but no new integrable cases were found. Probably only these three 
integrable cases exist among "diagonal" Hamiltonians. 

In the paper j^j the first integrable "non-diagonal" Hamiltonian (|1.2J) with k = (i.e. a Hamiltonian 
of the Kirchhoff type describing the motion of a rigid body in ideal fluid) was found. This Hamiltonian 
has a fourth degree additional integral. A generalization of this Hamiltonian to the case k ^ was 
reported in jHj. In the paper [S] the Hamiltonian has been rewritten in the form 

H = (M, AM) + (b, M x T), (1.3) 

where A is a constant symmetric matrix, b ^ is a constant vector and x stands for the skew product. 
It turns out that this class is very rich in integrable cases. In the paper [9 all Hamiltonians (|1.3|) with 
a quartic additional integral were described. Moreover, it was mentioned in piJJ EE E] that the general 
Sklyanin brackets |E] for the AAA-magnetic model lead to some integrable Hamiltonians of the same 
kind. 

The goal of our paper is a systematic investigation of integrable real Hamiltonians (jl.3|) and their 
inhomogeneous generalizations 

H = (M, AM) + (b, M x r) + (k, M) + (n, T), (1.4) 
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where k and n are constant vectors. 

It is remarkable that all known Hamiltonians ()1.3j) possessing additional polynomial integrals have 
also some linear partial integrals 



where u and v are constant vectors. Namely, for some constant vectors p and q the following relation 



holds. This implies that the corresponding equations of motion preserve the constraint P = 0. These 
linear partial integrals turn out to be factors of the polynomial integrals for all homogeneous Hamil- 
tonians considered in this paper. 

Because of this reason we start our study with a subsection 2.1 devoted to linear partial integrals. 
An interesting subclass of one-parametric families of Hamiltonians (jl.3|) arises there. In the next 
subsection 2.2 we apply the Kowalewski-Lyapunov test, which is well-known in the Painleve analysis, 
to find all possibly integrable families from this subclass. Because of a continuous parameter in the 
Hamiltonian this test becomes extremely efficient. The main result of this consideration is a new 
integrable Hamiltonian on so(3, 1) with an additional sixth degree integral. 

In section 3 we are dealing with inhomogeneous Hamiltonians (jl.4[) . As we claim in section 4, there 
are no Hamiltonians of the form ()1.3j) having additional integrals of degrees from 1 to 8 other than 
examples described in section 2. We have also verified that we found in section 3 all Hamiltonians 
(jl.4j) having additional integrals of degrees from 1 to 6. 

All computations for the paper have been done by the specialized computer algebra package Crack. 
It is designed to solve overdetermined polynomial differential and algebraic systems with an emphasis 
on extremely large problems. A major concern for any operations performed by the program is the 
complexity of resulting expressions. In ^3] an overview of the package and examples for its use in the 
classification of integrable systems are given. 

In this paper we present integrable Hamiltonians and the corresponding additional integrals in a 
general vector form, which is invariant with respect to the orthogonal transformations. For computa- 
tions with vectorial expressions extra code was written. 

2 Homogeneous integrable cases 
2.1 Linear partial integrals 

In this section we describe all Hamiltonians of the form (jl.3|) having linear partial integrals ()1.5|) . Re- 
lation (|1.6|) is equivalent to a system of bi-linear algebraic equations for coefficients of the Hamiltonian 
and components of vectors u, v, p and q. Below we present the result of our investigation of this 
system. 

One can check that there are two different possibilities, either case 1: v = q = 0, or case 2: vector 
q is equal to vector b from formula ()1.3|) . 






(1.6) 
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Case 1: Calculations show that any Hamiltonian having a linear partial integral with v = q = 
belongs to a class of "vectorial" Hamiltonians of the form 

H = a (a, b) |M| 2 + c 2 (a, m) (b, m) + (b,Mxr), (2.7) 

where b and a are constant vectors, q are constant scalars. In this case u = b and p = c 2 a x b. In 
other words, for Hamiltonian ()2.7|) we have 

{H, (b, M)} = c 2 (a x b, M) • (b, M). (2.8) 

In the next subsection we will present four different integrable Hamiltonians of the form ()2.7|) that 
possess additional integrals of degrees 1, 3, 4 and 6. 



Case 2: In the case q = b the following conditions 

(b, v) = (b, u) = 0, p = £b, 
where £ is a scalar, have to be fulfilled. 

Case 2a: If the vectors v and u are not parallel, then without loss of generality we may assume that 
b = u x v. It turns out that in this case £ = and u and v are arbitrary vectors such that (u, v) = 0. 
The Hamiltonian is given by 

H = i|u| 2 |M| 2 + - (u, m) 2 - ! (v, M) 2 + (u x v, M x r) . (2.9) 

The partial integral P = (u, M) + (v, T) satisfies the relation 

{H, P} = (u x v, r) • P. 
This Hamiltonian has the following additional integral of fourth degree: 

I = ^P|M| 2 — 2(v, M)(M, T)j • P. (2.10) 
This integrable case was found in a different non- vector form in [HUH]. A Lax operator is presented in 

m 

Case 2b: The other possibility is that the vectors v and u are parallel: 

(b, v) = 0, p = £b, u = r]v. 

It turns out that in this case £ = rj and k = rf. We see that a real linear integral exists only for the 
so(4)-version k > of bracket (jl.lj) . The Hamiltonian is given by the formula 

H = -2r](v, M^j (z, M^j + ^b,Mxrj, b = vxz (2.11) 
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where v and z are arbitrary constant vectors. The partial integral P = (v, 77M + T) satisfies 

{h, p} = (b, r]M + r) ■ p. 

This Hamiltonian has the following additional integral of fourth degree (cf. [T2*l HUj): 

I = (z, (rjM - r)|M| 2 + 2M (M, r)) ■ P. (2.12) 
The eigenvalues aij of the matrix A from (jl.2|) satisfy the following relations 

a 3 = 0, ai« 2 = — r / 2 |b| 2 . 
Diagonalizing the matrix A, we get the possible canonical form of the Hamiltonian ()2.11j) 

H = r]{cMl - -Ml) + M l72 - M 27l . 
c 

Although it is real only for the so(4)-bracket, after the renormalization crj — > c of the arbitrary 
parameter c we get the Hamiltonian (see 0) 

H = cM\ - -M\ + M l72 - M 27l , 
c 

which is real for any bracket In particular, if k = 0, we have a new integrable case on e(3) with 

a fourth degree integral. 

2.2 Kowalewski-Lyapunov test and a new integrable case. 

In this section we show that the class of Hamiltonians (|2.7jl contains a number of integrable cases. The 
first example of this kind was found in 

Example 1. Consider the so(3, Inversion k < of bracket (jl.lj) . The Hamiltonian 

= (a, b) |M| 2 - (a, m) (b, m) + (b,Mxr), (2.13) 

where the vector b is arbitrary and the length of the vector a = (01,02,03) is related to the Poisson 
bracket parameter k by 

a\ + a\ + a\ = -k, (2.14) 

possesses the additional quartic integral 



2 a, M x T - a, M - k\M 



b, M 



Recently in the paper the following integrable case has been found: 
Example 2. The Hamiltonian 



H = (a, b) |M| 2 - 2 (a, m) ^b, m) + (b, M x t\, (2.15) 



has under condition (J2.14j) the additional cubic integral 



b, MJ ^2 (a, M x Tj — k\M\ 2 + |r| 

In the two examples mentioned above, constraint (|2.14|) is necessary for integrability. We are going 
to find all integrable Hamiltonians similar to (|2.13|) and (J2.15|) . To do that we apply the Kowalewski- 
Lyapunov test to the class of Hamiltonians ()2.7|) assuming that the additional condition (j2.14|) is 
valid. 

Suppose we have a dynamical system 

^ = F(X), X={x u ...,x N ), F = (fx, . . . f N ) (2.16) 
at 

where are homogeneous quadratic polynomials of X. Solutions of the form 

X = j (2.17) 

for system (|2.1fiJ) with K being a constant vector are called Kowalewski solutions. Substituting 1)2.17)1 
into (|2.16|) . one obtains a system of algebraic equations for possible vectors K. 

The linearization X = Xo + of system (J2.16)) on a Kowalewski solution Xo obeys 

f = t s <*>- < 2 - 18 > 

where S is a constant iV x A^-matrix depending on the Kowalewski solution. 

Solutions of ()2.18|) have the form ^ = s t~ k , where k is an eigenvalue and s is an eigenvector of the 
matrix S. The number 1 — A: is called Kowalewski exponent. 

According to the Kowalewski-Lyapunov test, system ()2.16)) is "integrable" if for any Kowalewski 
solution all corresponding Kowalewski exponents belong to an a priori fixed number set A. The struc- 
ture of A is closely related to analytic properties of general solution for ()2.16j) . The usual choice A = 7h 
is associated with the requirement that the general solution should be single-valued. The latter is a 
standard assumption for the Painleve analysis. The most general version A = Q is associated with 
general solutions having algebraic branch points. But the main property for us is that A cannot be 
too wide. In particular, it cannot contain any open subset of C or R. Therefore for any one-parameter 
family of integrable (in the Kowalewski-Lyapunov sense) dynamical systems ()2.16|) the Kowalewski 
exponents must not depend on the parameter. This gives us strong necessary integrability conditions 
for one-parametric families of homogeneous quadratic dynamical systems. 

Given C\ and C2, the Hamiltonian ()2.7j) with (|2.14|) depends on one essential continuous parameter. 
Indeed, the length of a is fixed by (|2.14|) . the length of b can be normalized by the scaling of the 
Hamiltonian. The transformations 

M -> TM, r — > TV (2.19) 
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for any constant orthogonal matrix T preserve brackets (Jl.l)) and the form of Hamiltonian (J2.7j) . Two 
vectors of a fixed length have only one invariant (the angle between them) with respect to orthogonal 
transformation (|2.19|) . 

We want to find all pairs c\ and C2 in (|2.7j) such that the Hamiltonian is integrable for any value 
of the angle. 

Using transformation (J2.19|) . we may reduce b and a to 

b = (0,0,1), a=(a 1 ,0,a 3 ). (2.20) 

Taking into account the constraint (J2.14|) . we see that now a 3 remains to be the only free parameter 
in (J2.7)) . For a generic Hamiltonian of this kind the Kowalewski exponents for the equations of motion 
depend continuously on a 3 . For an integrable Hamiltonian these exponents must not depend on the 
parameter at all. 

Theorem 1. Suppose all Kowalewski exponents for Hamiltonian (|2.7|) with ()2.14|) . ()2.20|) do not 

depend on a 3 ; then the pair of constants ci, c% belongs (up to the transformation c\ — » — c±, c 2 — > — C2, 
which corresponds to a — > —a) to the following list 



a) 


Cl 


- arbitrary, 


c 2 = 


0; 


b) 


Cl 


= 1, 


c 2 = 


-2; 


c) 


Cl 


= 1, 


c 2 = 


-1; 


d) 


Cl 


= 1, 


c 2 = 


l . 

2' 


e) 


Cl 


= 1, 


c 2 = 


1. 



To prove this statement, we have to investigate the Kowalewski exponents on all solutions of the form 

Mi = -j-, 7* = f (2.21) 

for the equations of motion defined by Hamiltonian (J2.7|) . (J2.14|) . To get the list of Theorem 1, it turns 
out to be enough to investigate the Kowalewski exponents on special solutions ()2.21|) . which satisfy 
<73 = 1. There exist two classes of such solutions. For the first class we have 02(0.1^2 — 02^1) = 1- 
The second class is defined by m 3 = 0. For solutions of the first class, besides the cases of Theorem 1, 
the only extra surviving possibility is c\ — 1/2, c 2 = — 1. But this case does not pass the Kowalewski- 
Lyapunov test on the solutions of the second class. For each pair (c\, 02) of the remaining list a)-e) of 
Theorem 1 we verify that for all Kowalewski solutions the Kowalewski exponents do not depend on 
a 3 . All computations are straightforward. A technical problem is that the explicit form of solutions of 
the first class involves radicals. To avoid this difficulty, we used calculations based on the Groebner 
basis technique. 

Comments. The Hamiltonian a) belongs to the family 

H = ci|b| 2 |M| 2 + c 2 (b, m) 2 + (b, M x t), (2.22) 
which possesses (without restriction (|2.14j0 the linear integral of motion / = (b, M). 
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The Hamiltonians b) and c) are described in Examples 2 and 1, correspondingly. The possibility 
d) leads to the following new integrable case with an additional integral of sixth degree: 

Example 3. The Hamiltonian 

H = (a, b) |M| 2 - - (a, m) (b, m) + (b,Mxr), (2.23) 
under condition (|2.14jl has the additional sixth degree integral 




-r 2 (M,bxa) -(bxa,Mxr) -M 2 r 2 (b,a) - KM 2 T 2 b 2 . 

Possibly this integrable case is related to the simple Lie algebra g{2). 

In [T3| a generalization of the general scheme by Sklyanin ^2] has been proposed. In JT] the 
authors found a separation of variables for Hamiltonian (j2.15|) from Example 2 in the framework of 
this approach. Probably a separation of variables for the Hamiltonians from Examples 1 and 3 could 
be found after some development of these ideas. 

All Kowalewski exponents for cases a)-c) are integers. In case d) there are two solutions of the 
class 2. For these solutions we have 

det(S-fcld) = (k- l) 3 {k + 2){k- -)(k + |) 

and 

det(5 - k Id) = (k - 1) 3 (A; + 2){k + i) 2 . 

In other words, some of the Kowalewski exponents are half- integers. 

Case e) is a mysterious one. We have verified that the Hamiltonian has no polynomial additional 
integrals of degrees not greater than 8. On the other hand, on all Kowalewski solutions all Kowalewski 
exponents are integers. It would be interesting to verify whether the equations of motion in the case 
e) satisfy the standard Painleve test. 

3 Inhomogeneous integrable Hamiltonians. 

3.1 Admissible linear terms for integrable homogeneous Hamiltonians 

In this subsection we find for integrable homogeneous Hamiltonians H of the form (|2.22|) . (j2.15|) . (j2.13D . 
(I2~TT1) . (P or <m possible linear terms 

T=(k, M) + (n, T), (3.24) 
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where k and n are constant vectors, such that the Hamiltonian 

H = H + T 

has an additional integral of the same degree as H. 

Proposition 1. The following linear terms are admissible (in the above sense): 



1) 


for 


Hamiltonian (jZZD 


T 


= Pl (h, M)+p 2 (b, r); 


2) 


for 


Hflmiltrmifln (19 1 with ID 1 dt\ 


T 




3) 


for 


Hamiltonian fl37TT5fc with (HH1) 


T 


= (pia + Mxb, M)+p 3 (b, T 


4) 


for 


Hamiltonian (j2.11j) 


T 


= (k, M)+ Pl (vxz, r); 


5) 


for 


Hamiltonian ()2.9|) with (u, v) = 


T 


= Pi(u, M) + p 2 {u x v, r); 


6) 


for 


Hamiltonian (jZ2SD with (|2~T2|) 


T 


= pi(a x b, M), 


where k 


is an arbitrary vector and Pi,P2, P3 


are 


arbitrary constants. 



We present the explicit form of the additional integrals for the non-homogeneous Hamiltonians H 
of Proposition 1 in appendix A. 

3.2 A deformation of the Poincare model 

The formula ()2.22j) describes Hamiltonians (jl.3j) that have linear additional integrals. We shall call 
()2.22|) the Poincare model. 

For the special pair c\ = 1 , c 2 = — 1/2 the Poincare model is superintegrable. It means that besides 
the linear integral the Hamiltonian has a polynomial integral of degree higher than 1. 

Proposition 2. The Hamiltonian 

H hom = |b| 2 |M| 2 - i(b, M) 2 + (b, M x r) (3.25) 

under condition |b| 2 = — k has the following additional integral of degree 4, functionally independent 
of H, the Casimirs and the linear integral (b, M): 

/ = (k, M) [(k, M)|b| 2 — 2(k, b)(b, M)] ■ [|T| 2 + |b| 2 |M| 2 ] 
+ |M| 2 (b, M) 2 [(k, b) 2 + |k| 2 |b| 2 ] - (kxb.Mxr) 2 
+2^k x b, b x m) ■ [|M| 2 (k,b x T) - (M, T)(k,b x M)] , 
where k is an arbitrary constant vector. 
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It follows from the Jacobi identity that I\ = {(b, M), 1} is a first integral as well. It turns out that 
I\ 7^ 0, which means that the Poisson subalgebra of polynomial integrals for ifhom is non-commutative. 
For any X let us denote the Poisson bracket {(b,M), X} by X'. One can check that the integral / 
satisfies the relation I'" + 4|b| 2 /' = 0. 

The derivation X — > X' can be regarded as a linear operator on the finite dimensional vector space 
S n of all homogeneous polynomials of degree n depending on components of M and F. The operator 
spectrum is /i^ = zA;|b|, 1 < k < n. Probably this operator plays an important role in the theory of the 
vector Hamiltonians ([2.70 . Many terms in the integrable Hamiltonians from this class admit simple 
descriptions in terms of this operator. For example, the linear polynomial (k, b x M), where k is an 
arbitrary constant vector, is the general linear solution of the equation X" + |b| 2 X = 0. This fact and 
formula ([2.8(1 imply that the arbitrary Hamiltonian (j2.7[) satisfies the equation X'" + |b| 2 X' = 0. The 
additional integral from Example 1 satisfies the same equation and so on. 

The Hamiltonian ([3.25(1 admits the following inhomogeneous integrable extension: 

H = H hom + (kxh, M)+pi(b, T), (3.26) 

where k is an arbitrary constant vector, pi is an arbitrary constant. Hamiltonian ([3.26(1 under condition 
([2.14(1 has an additional integral of degree 4, given in the appendix A. 

4 Classification results 

It is very likely that all integrable Hamiltonians of the form ([1.3(1 and ([1.4(1 are exhausted by the 
examples presented in sections 2 and 3. 

Theorem 2. Suppose a Hamiltonian of the form ([1.3(1 with real coefficients has an additional 
polynomial integral of degree from 1 to 8; then the Hamiltonian belongs to the six families ([2.22(1 . 
(|2~T5) . (jUSD , (Hill), O or (jUj . 

Scheme of the proof. Using transformations ([2.19(1 . one can reduce any (real) Hamiltonian ([1.3(1 

to 

H = a x M\ + a 2 M 2 + a 3 M 3 2 + a 4 MiM 3 + a 5 M 2 M 3 + M^ 2 - Af 2 7i. (4.27) 

In this canonical form the vector b is normalized to (0,0,1). The alternative idea of bringing the 
matrix A to the diagonal form lead to overwhelming computational complexity. 

Without loss of generality we may assume that the additional polynomial integral is homogeneous. 
Given the degree m of the additional integral /, we form the general homogeneous m-th degree poly- 
nomial of six variables M$, 7$ with undetermined coefficients. The condition {I, H} = gives rise to 
a bi-linear system of algebraic equations for both coefficients of H and I. Of course, this system can 
be solved "by hand" only for small m. If m = 6 the algebraic system contains 791 bi-linear equations 
for 458 unknown coefficients. This system can not be currently solved by standard computer algebra 
systems. Also all attempts by the authors to use the two best known packages specialized in the 
solution of polynomial algebraic systems failed. 
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The computation was performed using the computer algebra package Crack. For m < 6 the 
calculations were performed automatically and for m > 6 with manual interaction. 

Following the same line we have obtained 

Theorem 3. Suppose a Hamiltonian of the form (jl.4j) with real coefficients has an additional 
polynomial integral of degree from 1 to 6; then the Hamiltonian belongs to the seven families described 
in Propositions 1 and 2. 

5 Conclusion: unsolved problems 

This paper as well as [5] belongs mostly to so called "experimental" mathematical physics. The result 
of the experiment is a new interesting class (jl.4|) of quadratic Hamiltonians. This class contains 
several new integrable cases. Theorems 1-3 give reasons to believe that we found all real integrable 
Hamiltonians of the form ([1.4)1 . However, this should be proved. To do that, one can apply the 
Painleve approach or methods developed in [To] . 

The separation of variables for several models from our paper is also an open problem. 

Besides Hamiltonians listed in Theorem 3, there exist integrable Hamiltonians of the form (|1.4j) with 
complex coefficients. For example, the Hamiltonian (J2.7j) . (J2.14)) with ci = 1, c-i — — | has an additional 
integral of sixth degree under condition |b x a| = 0. To find the complex Hamiltonians, one should 
consider two different normalizations of the vector b. The first is b = (0,0,63), which corresponds 
to the possibility |b| 7^ and gives rise to the normal form ([4.27)1 . The second normalization b = 
(0, 62, i 62) corresponds to |b| = 0. In the case of fourth degree additional integrals all complex integrable 
Hamiltonians have been found in [Oj. 

Probably all Hamiltonians from our paper have their quantum counterparts (for the case of fourth 
degree integral see [9 ). It would be interesting to find the corresponding quantum Hamiltonians and 
integrals as well as the quantum separation of variables. 
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Appendix A 

Here we present an explicit form of additional integrals for inhomogeneous Hamiltonians from Propo- 
sitions 1 and 2. Each integral I of degree m is a sum / = Y^Li where U is a polynomial of degree i 
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homogeneous in M, T. 

Case 1 of Proposition 1. I = I\ = (b,M). 

Case 2 of Proposition 1. I 3 is presented in Example 2, 



h = 2 Pl (b,M)(a,r)-(k,a)|M| 2 -(k,Mxr), 

h = -#,M)- Pl (k,r), 



Case 3 of Proposition 1. J 4 is presented in Example 1, 



2(b,M) \ Pl (a,M) 2 + /t|Mp 



+ P2 



2(a,M x T) 



(b x a, M x r) + (a,M)(a,b x M) + p 3 (b, M)(a, T)} 



2|r| 2 - (a,M) 2 + 2(a,M x V) - 2p 2 p 3 (b, M)(a, b x T) 

b| 2 (a, M) 2 - «(b, M) 2 - 2(a, b)(a, M)(b, M)j - 4pip 3 (b, M)(a, T) 
p 2 3 K(b, M) 2 + 2p x p 2 [(a x b, M x T) - (a, M)(a, b x M) 



+pl 



h = 2p x p 3 



Pa (a, T) + p 2 (a, b x T) + p 3 /c(b, M) 



Case 4 of Proposition 1. J 4 = 77 1 a x b| 2 / where I is given by (|2.12jl . 



pi^l (v, z) 2 - |v| 2 |z| 2 



77|M| 2 (v, z x r) - 2(v, r)(z, m x r) + |r| 2 (v, z x m) 



+r] 2 \M\ 2 (z, M)(v x z, k x v) - (v, M)(v x z, k x z) 



+2r] 



(v, M)(z, k)(v x z, T x M) + (v, M)(k, M)(v x z, z x T) 



-(v, M)(k, r)(v x z, z x M) 

~(k, r)|M| 2 |v x z| 2 - (z x M, M x T)(v x z, v x k) 
-(v,z x k)(v,z x M)|r| 2 , 



p\r]\w x z| 2 rf(y x M, z x M) + (v, T)(z, V) - p x r) ^(v, M)(z, V) 
+ (v, T)(z, M)) (v, z x k) + (v x z, z x k)(v, MxT) 

+ (v x z, v x k)(z, M x r)l - ??|M| 2 (v x k, z x k) + (v, z x k)(k, M x T), 



h = Pl (v, z x k) • (W(v, z x M) + (k, r)) . 
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Case 5 of Proposition 1. 7 4 is is given by (|2.1(Jj) . Ii = 0, 

J 3 = 2 Pl (p|M| 2 -(v,M)(M,r))+2p 2 ((M,r)(u,vxM)-P(M,rxv)), 

h = ^|M| 2 +p^| v |2|r|2 + fi . (V)M) 2_ (V)r)2 j _ 2m ( M) rx v ), 

where P = (u, M) + (v, T). 

Case 6 of Proposition 1. J 6 is presented in Example 3, l\ = 0, 

h = 4pi(b,M) | (a,b)|M| 2 (a x b, a x M) (a, b x M) + |T| 2 (a x b, b x M) (a, b x M) 
+ (a x b, M x T) [(b, T)(a, b x M) - (b, M)(a, bxT 
+ (b, M)(M, T) (a, b)(a x b, a x M) - «(a x b, b x M) 
-(a, b)(a x b, b x M) |2(a, M)(M, T) - (a, r)|M| 2 
+ |M| 2 (a x b, a x M) |b| 2 (a, T) - 2(a, b)(b, T) } , 

h = pj {4 (a x b, b x M) ((a, T)(b, M) + 2 (a, b)(M, r)) - (b, T)(b, M) (a x b, a x M) 
-2(a, b)|M| 2 (a x b, b x T)j • (a, b x M) + 4(b, M) [(b, M)(a x b, a x M) 
-(a, M)(a x b, b x M) ■ (a, b x T) - 8(a x b, b x M)(a x b, b x r)(M, T) 
+ 4(a, b) 2 (a, M)(a x b, b x M) - 4(b, M) ^2(a, b) 2 - |a| 2 |b| 2 ) (a x b, a x M) 
+ 4(a x b, b x T) 2 ] |M| 2 - 4(a, b) 2 ((a, b) 2 - |a| 2 |b| 2 ) |M| 4 } , 

J 3 = 8p?|a x b| 2 (^(a,b)|M| 2 (a,b x M) + (b x a, b x T)|M| 2 + (a x b, b x M)(M,r 



J 2 = -4^|ax b| 2 ((a,b) 2 |M| 2 - |b| 2 |r 12 



The Case of Proposition 2. Ii = 0; J4 is related to integral / from Proposition 2 by 
h = |b| 2 J + (b,M) 2 [(k x b, b x T) (k,b x M) - (k x b, b x M) (k,b x T) 



+ ((k, b) 2 -|k| 2 |b 12 



1 



-(b,M) 2 + |k| 2 |b| 2 [\T\ 2 + k\M\ 



Pi 



({k, b) 2 - |k| 2 |b| 2 ) (b, M) 2 (b, r) - 2|b| 2 (kxb.Mxf) ((k, b x T) + (k x b, b x M)) 



+ (k, b) 2 - |k| 2 |b| 2 (b,M) (b,M) (k,b x M) + 2(k x b, T x M) 
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h = Pi |b| 2 (k x b, b x M) 2 + (k x b, b x T) 2 + ((k, b) 2 - |k| 2 |b| 2 )|b x if 
+ |k| 2 |b| 4 (|r| 2 + /t|M| 2 )] - 2 Pl ((k, b) 2 - |k| 2 |b| 2 ) (b, M) (k, b x r) + 
((k, b) 2 - |k| 2 |b| 2 ) (b, M) (k x b, k x M) - (k, M) (k x b, b x M)] . 

Appendix B 

With the comments in this appendix we only want to give some impression on features of the package 
Crack that we used for the solution of the overdetermined algebraic systems {/, H} = 0. The 
description is neither complete not detailed. 

1. Flexibility (different possibilities to change or adapt the strategy): 

a) CRACK has a priority list, which contains among other things, 10 different types of substitu- 

tions, two types of factorizations, three types of shortening of equations, various types of 
Groebner basis calculation steps and the possibility to call external specialized packages. 
One can set the different priorities for these internal procedures before the session but also 
change them interactively during the computation. In addition, the program is able to 
make minor changes in the priority list in accordance to its analysis of the current situation 
and history of the session. An important point for the analysis is that each equation has a 
property list describing the previous use of the equation as well as its potential for further 
computation. 

The program has automatic and interactive modes and the possibility to switch between 
both during the computation. In the interactive mode it provides extended support for 
inspecting the system of equations, to specify the very next steps in detail or modify the 
solution strategy in general. 

b) After a decision is made what kind of step should come next (based on the priority list) more 

heuristic knowledge decides how to do it. For example, if a factorization is to come next 
then a refined weighting scheme determines which equation to factorize and in which order 
to set its factors to zero. 

c) During computation inequalities are actively gathered and simplified in order to reduce the 

number of branches. 

2. The program is paranoid about expression swell and for its prevention it is able to 

a) compute a tight upper bound for the system's length after an intended substitution and hence 

choose substitutions giving minimal growth, 

b) find linear combinations of pairs of equations in an attempt to shorten and simplify the 

system recursively during computation, and 
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c) it takes advantage of the initial bi-linear form of the system and keeps it linear in the unknown 
coefficients of the first integral at all times. 

3. Several facilities are provided, e.g. 

a) to find re-parametrizations of solutions in order to merge them, 

b) safety precautions: to interrupt intermediate steps automatically if they (unexpectedly) take 

too much time, to catch all interactive input, to conveniently store and load backups, 

c) an automatic preparation of web pages with results of the computation. 

4. Different branches of the general solution can be investigated at the same time in parallel on a 

cluster of computers. 

To obtain a free copy of this REDUCE package contact T. Wolf under twolf Obrocku. ca . 
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